Statistical Anisotropy and the 
Vector Curvaton Paradigm 

Konstantinos Dimopoulos^ 

Physics Department, Lancaster University, Lancaster LAI 4^B, UK 

,-H Abstract 
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^SJ ■ The vector curvaton paradigm is reviewed. The mechanism allows a massive vec- 

^ . tor boson field to contribute to or even generate the curvature perturbation in the 

Universe. Contribution of vector bosons is likely to generate statistical anisotropy 
in the spectrum and bispectrum of the curvature perturbation, which will soon be 

^^ . probed observationally. Two specific models for the generation of superhorizon spec- 

tra for the components of an Abelian vector field are analysed. Emphasis is put 
on the observational signatures of the models when the vector fields play the role 
of vector curvatons. If future observations support the vector curvaton mechanism 
this will open a window into the gauge field content of theories beyond the standard 
model. 
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(N ■ 1 Introduction 
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S: 



Cosmic inflation is arguably the most compelling way to overcome, or at least ameliorate. 



C"^ \ the so-called horizon and flatness problems of the hot big bang cosmology. However, these 

. ■ problems are successfully addressed by all models of inflation, provided that the inflationary 

Q ■ expansion lasts long enough. Therefore, discrimination between inflation models is based 

on another important consequence of inflationary expansion, namely the generation of 
the curvature perturbation C, in the Universe, which is responsible for the formation of 
structures such as galaxies and galactic clusters and which is reflected onto the Cosmic 

^ ■ Microwave Background (CMB) radiation through the Sachs- Wolfe effect. 

The latest CMB observations appear to conflrm the "vanilla" predictions of inflation 
with respect to C,: scale-invariance, Gaussianity and statistical homogeneity and isotropy. 
However, a period of accelerated expansion of space (this is the deflnition of cosmic infla- 
tion) is not enough to guarantee scale-invariance for the curvature perturbation. Indeed, 
inflation is required to be of the quasi-de Sitter type, where the density of the Universe 
remains roughly constant. Also, the CMB observations seem to suggest that exact scale- 
invariance is not favoured (although it is not ruled out) with the spectral index of C? 
satisfying n^ — 1 = —0.037 ± 0.014 (at la) [1] when ACDM cosmology is assumed.^ Most 
likely, this reveals some dynamics for inflation (the density is not exactly constant), which 
is indeed expected by model-builders. Similarly, a high degree of Gaussianity in the curva- 
ture perturbation is expected since this reflects the randomness of quantum fluctuations, 
on which the particle production process is feeding, in order to generate the superhorizon 
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^Exact scale invariance corresponds to ti^ = 1. 



spectrum of the fields which eventually give rise to (. However, the Gaussianity of ( cru- 
cially depends on the linearity of the process which translates these field perturbations to (. 
This is quantified through the so-called non-linearity parameter /nl- The latest CMB ob- 
servations provide a hint of non-zero non-Gaussianity since, in the squeezed configuration, 
/nl = 32 ± 21 (at la) [1], which again appears to deviate from the "vanilla" prediction.^ 

In the same spirit, observations suggest that there may be deviations from statistical 
homogeneity, since there seem to be a difference in the power of ( as large as 10% between 
hemispheres in the CMB [2]. Finally, statistical isotropy in ( is also questioned by obser- 
vations. Indeed, there is tantalising evidence of a preferred direction on the microwave sky. 
This is the so-called "Axis of Evil" observation [3], which amounts to an alignment of the 
quadrupole and octupole moments in the CMB, which is statistically extremely unlikely 
[4] and has been shown to persist beyond foreground removal [5]. 

Thus, we see that the precision of the cosmological observations is such that begins 
to enable us to explore beyond the "vanilla" predictions of inflation and use the observed 
deviations from them to discriminate between classes of models and paradigms. There 
is already a huge literature on deviations from exact scale-invariance and Gaussianity for 
the curvature perturbation. In this paper we discuss possible deviations from statistical 
isotropy and we present a compelling paradigm for their generation; namely the Vector 
Curvaton Paradigm. 

Throughout the paper we consider a metric with negative signature and use natu- 
ral units where c = h = ks = 1 and Newton's gravitational constant is SttG = rrip^, with 
nip = 2.4 X 10^* GeV being the reduced Planck mass. 

2 Statistical anisotropy in the curvature perturbation 

Statistical anisotropy amounts to direction dependent patterns in the CMB. It can be 
quantified as follows. The spectrum Vq of the curvature perturbation is defined through 
the two-point correlator as 



where k = \k\ and 



C(fc) = fc{x)e-'''-''d^x. 



The reality condition C*{k) = C{~k) demands that V({k) = V(^{—k). Now, the dependence 
of the power spectrum on the direction of the momentum vector can be parametrised as [6] 

V^{k)=Vf{k)[l+g{d.kf + ...l (2) 



■^Note that Q is indeed predominantly Gaussian because the bispectrum Bq is related to the power 
spectrum Pc^ roughly as Bc_ ^ JnlP^ (see Eq. (15)), with the observations suggesting Vq = 4.8 x 10~^ for 
the spectrum [1]. 



where d is the unit vector along the preferred direction, k = k/k and the elhpsis denotes 
higher than quadratic order terms, which are neghgible ii g < 1. g = g{k) is sometimes 
called the anisotropy parameter. A similar parametrisation can be assigned to higher order 
correlators, i.e. the bispectrum, trispectrum etc. (for the bispectrum see Sec. 3.2). 

What are observations saying about g7 In Ref. [7] the latest CMB data was analysed 
without the prior of statistical isotropy, in order to obtain bounds or even detect statistical 
anisotropy. Indeed, it was found that g = 0.29 ± 0.03 at the level of 9a\. However, the pre- 
ferred direction was too close to the ecliptic plane so the authors suspected some unknown 
systematic.^ Hence, this number can be considered only as an upper bound g < 0.3. The 
observations of the Planck satellite will strengthen this bound by at least an order of mag- 
nitude and reduce it to (? < 0.02 if statistical anisotropy in the spectrum of the curvature 
perturbation is indeed not observed [9]. 

3 Vector fields and the curvature perturbation 

Statistical anisotropy in the curvature perturbation cannot be generated if one considers 
the effects of scalar fields only, because the latter cannot introduce a preferred direction. 
In this paper we will study how vector boson fields can directly influence the curvature 
perturbation and generate statistical anisotropy.^ We will investigate the contribution of 
vector fields to ( through the so-called 6N formalism [11]. 

According the 5N formalism the curvature perturbation is the difference of the logarith- 
mic expansion between uniform density and spatially flat slices of spacetime: ( = 6{ln a) = 6N, 
where a is the scale factor of the Universe and A^ corresponds to the elapsing e-folds of 
expansion. We will assume that A^ is influenced by both scalar and vector boson flelds. For 
simplicity, we consider one of each of such flelds^, i.e. A^ = N{(f), A). Then, the curvature 
perturbation can be written as 

C = Ar^50 + N\6A, + ^N^^{5<PY + ^NiJ<pSA, + ^N'/JAMj + ■ ■ ■ , (3) 

where N, ^%,N\^ |£, N,, = AT;^ = ^ and Nf^ ^ ^, with z = 1,2,3 
labelling spatial components and Einstein summation over repeated indexes is assumed. 
Now, since the vector fleld has three degrees of freedom, at a flat slice of spacetime foliation, 
we deflne 

5A{k,t) = Y.e^{k)5A^{k,t), (4) 

A 



"'it is important to note that this finding was also confirmed by the WMAP team in Ref. [8], who 
suspect that it is due to beam assymmetries but have not fully expalined it yet. 

^Indirectly, statistical anisotropy in C, can also be generated by considering a mild anisotropisation of 
the inflationary expansion, due to the presence of a vector boson field condensate. In this case, it is the 
perturbations of the inflaton scalar field which arc rendered statistically anisotropic [10]. 

^For the multifield case see Ref. [12]. 



where X = L,R, \\ denotes the three polarisations and the polarisation vectors can be de- 
fined as 

eL = ^(1,1,0), en = ^{1,-1,0) and e|| = (o,0,l), (5) 

where 'L', 'R' denote the left and right transverse polarisations respectively and '||' de- 
notes the longitudinal polarisation (if physical). Then, assuming approximately isotropic 
expansion^, the power-spectrum for each polarisation of the vector field perturbations is 



27r2 



{6A^{k)6A^{k')) = {27rf6{k + k')—Vx{k) 



(6) 



3.1 The spectrum 

In Ref. [13], it was shown that the correlators of the perturbations of the vector field can 
be written as 



{6Ai{k)SAj{k')) = {27if6{k + k 



.2tt^ 



A;3 



7;t(fc)p+ + zT,7(fc)p_ + Ti(fcm 



(7) 



where 



7^+(fc) = 6ij - kikj , T. . (fc) = Eijkkk and TlUk) = kik 



hu^ 



(9) 



with 6ij being the Kronecker's delta, and we have defined for the transverse spectra 

denoting the parity even and odd polarisations. For a parity conserving theory we have 
V- = 0. Using the above, we obtain the power spectrum of the curvature perturbation as 

Nlvm + N\N^A \n{k)V4k) + Tl{k)V\\{k) = 



Vdk) 



NiVs + m 



(10) 



V+ + {V\\-V^){NA-kf\, 

where Na = \Na\ = ^ N\N\, Na = N^/Na. 

Prom the above we see that the isotropic part of the spectrum is 

V^\k) = NlV^{k) + NlV+{k) (11) 

and the preferred direction is given by cZ = Na (cf. Eq. (2)). The anisotropy parameter is 
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/3 






where we have defined 



Nl 



/3 = -^ 



(12) 



(13) 



which quantifies the relative contribution of the vector over the scalar field to the modu- 
lation of A^. Notice that particle production becomes isotropic {g = 0) if V+ = V\\. 



^This is in contrast to Ref. [10] where it is the anisotropy in the expansion which sources statistical 
anisotropy in C,. 



3.2 The bispectrum 



Statistical anisotropy is also possible to be manifest in higher order correlators of the 
curvature perturbation. In this paper we discuss only the bispectrum (for the trispectrum 
see Ref. [14]). 

The bispectrum of the curvature perturbation is defined as 



{ak)ak')ak")) = {27rf5{k +k' + k")B^{k, k\ k") . 



(14) 



The bispectrum Bq is a measure of the non-Gaussianity of the curvature perturbation since, 
for Gaussian C,, B^ is exactly zero. 

The curvature perturbation is generated due to the quantum fluctuations of suitable 
fields which are stretched to become classical perturbations during inflation. Since quantum 
fluctuations are Gaussian (which reflects their randomness) sizable non-Gaussianity in ( is 
generated only if the process through which the perturbations of the relevant fields affect 
the expansion of the Universe and imprint their contribution to the curvature perturbation. 
If this process is significantly non-linear deviations from Gaussianity will be generated. This 
is why the bispectrum is quantified by the so-called non-linearity parameter /nl, which 
can be defined as follows 

Bdki, k2, ks) = ^/nl [Pdh)P^{h) + P^{k2)Pc{h) + P^{h)Pc{h)] , (15) 

where iirk^PQ = {2txYVc_. The value of /nl depends on the configuration of the three mo- 
mentum vectors which are used to define the bispectrum. The most popular configurations 
are the "equilateral", for which ki = k2 = k^, and the "squeezed", for which ki = k2 ^ k^. 
How does the contribution of a vector field affect the bispectrum of the curvature 
perturbation? In Ref. [15] it was shown that 



Br 



Bd, + Baa + B/ 



where 

Ba 



nINa 



47r^ A-K^ 



471^ 
J,3l,3 



V4,{hi)V^{ki) 



(16) 



(17) 



B, 



(f>A 



-^n^nIa 



■pZa 

"'l"'2 



V^{ki)Aii{k2) + 5 cyclic permutations 



and 



5A=^A^.(fcl)A^ilA^,(fc2) + ^M(fc2)i^A-^.(fc3) + ^M(fc3)iVi4-^,(fcl), (19) 



U31.3- 



where 



H2H3 H3H1 

M{k) = V+{k)NA \Na + p{k)k{k ■ Na) + iq{k)k X Na 



(20) 



and we have defined 

V\\ - V+ V- 

P = — ^ and q = — — . (21) 

Using tlie above, we obtain /nl in the equilateral and squeezed configurations respectively 
as follows 

6 , ^ Bf{k,MM) 



and 



6 Bl-\k„k„k,) 



, 2 „„, , , 1,31,3 



where B^^ = (^) -B^ and B^^^ = ^^-B^*^^, with ki = k2^ k^ in the squeezed configu- 
ration and A; = fci = /c2 = ^3 in the equilateral configuration. 

4 The Vector Curvaton Paradigm 

For a vector field to directly affect the curvature perturbation in the Universe we need two 
ingredients. First, we need a mechanism to break the conformal invariance of the vector 
field and generate a superhorizon spectrum of vector field perturbations 6Af^. Second, 
we need a mechanism that will allow these perturbations to affect (or even generate) the 
curvature perturbation (. This can be done only if the vector field and/or its perturbations, 
in some way affect the Universe evolution. 

In this section we focus on the second ingredient, i.e. on a mechanism for the generation 
of a contribution from the vector field perturbations to the curvature perturbation of the 
Universe; namely the Vector Curvaton mechanism. Thus, we assume that some given 
mechanism has produced the necessary superhorizon spectrum of vector field perturbations 
(as is discussed in Sec. 5) during infiation, which for the moment we take for granted. 

A single vector field cannot play the role of the infiaton. The reason is straightforward. 
Infiation homogenises a vector field and a homogeneous vector field picks up a preferred 
direction in space. ^ Thus, if a homogeneous vector field dominated the Universe during 
infiation it would lead to excessive anisotropic stress, which would produce too much of a 
large-scale anisotropy and, therefore, will be in confiict with CMB observations. A huge 
number A/" of vector fields, randomly oriented, could avoid this problem [17]. Indeed, if 
this is the case then the statistical anisotropy produced is (7 oc I/VAT, which means that 
hundreds of vector fields are needed to satisfy the observational bounds. This not only 
implies the use of giant gauge groups but also requires the tuning of the initial conditions 
so that they are the same for all the fields. Another option is to consider a "triad" of 
orthogonally oriented vector fields (again with the same initial conditions) so that the 
excessive anisotropic stress is eliminated [18]. For the above reasons we will not consider 
vector fields as infiatons. 



^Unless one tunes the spatial components of the vector field to zero by design [16] 



If the vector field is not the inflaton, it needs to affect the Universe expansion in 
some other way, either at the end or after the end of inflation. There are a multitude 
of mechanisms which may allow a vector field to do that, mirroring the corresponding 
scalar field models. Prominent examples include the curvaton [19], the inhomogeneous 
end of inflation [20] (see also Ref. [21]) and the modulated reheating [22] mechanisms. 
Historically, statistical anisotropy by vector field perturbations was first studied in the 
context of the inhomogeneous end of inflation mechanism [23] ,^ using a particular model of 
hybrid inflation. ^° Here, however, we concentrate on the curvaton mechanism, which has 
the considerable advantage that it does not rely on an interaction of any kind between the 
vector field and the inflaton sectors. This avoids overcomplicating the model but, more 
importantly, allows the vector curvaton sector to be completely independent of the physics 
of inflation. As such, not only can it apply in many given inflation scenarios, but it can 
correspond to physics at a much lower energy scale than inflation; even TeV physics. As 
is the case of the scalar curvaton, the vector curvaton is not a particular model but it can 
correspond to a multitude of realisations, hence we refer to the mechanism as a paradigm 
rather than a model. The vector curvaton mechanism was first introduced in the pioneering 
work in Ref. [26] , which was the first article to consider the possibility that a vector boson 
field can contribute to the curvature perturbation in the Universe. 

The idea of the curvaton assumes the existence of a spectator field during inflation, 
which has nothing to do with inflationary dynamics but is light enough so that it manages 
to obtain a superhorizon spectrum of perturbations. After the end of inflation (possibly 
long afterwards), the curvaton becomes heavy and begins undergoing oscillations which 
allow it to come to dominate (or nearly dominate) the Universe before its decay. Owing 
to its perturbations, the density of the curvaton is perturbed throughout space so that its 
(near) domination occurs at different times at different locations. Thus, its effect on the 
evolution of the Universe is location dependent, which is the reason why it can affect (or 
even generate) the curvature perturbation in the Universe. Note that, for a vector field to 
do this, it must avoid generating an excessive anisotropic stress at domination. 

4.1 The setup 

Consider a massive Abelian vector boson field, with Lagrangian density 

£ = -^F,..F^" + ]^m'W^W^, (24) 

where W^ is the vector field, m is its mass and F^^, = d^W^ — di,W^ is the field strength 
tensor. Inflation homogenises the vector field so that W^j_ = W^{t). If ?7i 7^ it is easy to 
show that the temporal component of the homogeneous vector field is zero, i.e. Wt = 0. If 
m = then the field is gauge invariant and we can set Wt = hy virtue of a gauge choice. 
However, in this case the value of the spatial vector field W is not well defined because 



^For modulated reheating with vector fields see Ref. [25]. 
^°For non-Gaussianity in this model see also [24]. 



gauge invariance allows us to change it as W^W+C, where G is a constant vector of 
arbitrary magnitude. Thus, we will concentrate on the case m ^ from now on, where 
gauge invariance is broken and the homogeneous "zero-mode" is well defined. 
The energy-momentum tensor for the vector field is 

T^u = Ig^.uF.^F^'^ - F^pF^P + m^ (w,W, - ^g^.W.W^ , (25) 

where Qf^^, is the metric tensor (negative signature is assumed). The above can be written 

as [26] 

^/T = diag(pA, -P±, -Pu +Pl) , (26) 

where 

Pa = Pkin + Va and p± = pkin - Va (27) 

with 

PMn = -\F^.F^^'' and VA = ~m^W^W'', (28) 

where 'kin' denotes the kinetic density pkin- Notice that the energy-momentum tensor is 
similar to the one of a perfect fluid with the crucial difference that the pressure along the 
longitudinal direction is of opposite sign to the pressure along the transverse directions. 
This means that, if this pressure were not zero and the vector field dominated the Universe, 
it would give rise to significant anisotropic stress, which is the reason why a single vector 
field cannot play the role of the inflaton. 

Using Eq. (24) one can obtain the equation of motion for the homogeneous vector field, 
which reads 

W + HW + rn^W = , (29) 

where the dot denotes the derivative with respect to the cosmic time t and H = a/a is 
the Hubble parameter, i.e. the rate of the Universe expansion. At this point we need to 
stress that W is the comoving and not the physical vector field. Indeed, the mass term in 
Eq. (24) can be written as 

5Cm. = ^m^W^W = ^m^iW^ - a-^WiWi) = -Kn'\W/a\\ (30) 

where we used that VFj = and a spatially flat FRW metric ds^ = dt^ — a^dx^dx^. From 
the above it can deduced that the physical vector field has spatial components 

A = W/a. (31) 

In terms of the physical vector field, Eq. (29) is written as 

A + ?,HA + {H + 2H^ + m^)A = , (32) 

where A = |A|. From Eqs. (28) and (31) one finds 

Pi.in = liA + HAf and ^4 = ^mM^. (33) 

8 



The solution of Eq. (32) is of the form [27] 

A = t^(^) [ciJdimt) + C2 J„d(mt)] 



l+3w 



(34) 
ci,C2 are constants of 



where w is the barotropic parameter of the Universe, d _ .^^ , . 
integration and Jd denotes Bessel function of the the first kind with order d. 

When the physical vector field is light m '^ H ^^ mt ^ 1 the above solution can be 
approximated as [27] 



A 



f a 



2w + 1 Vacnd 



k3t«-l) 



^cnd + 



A. 



end 



H, 



end , 



(35) 



where the subscript 'end' denotes the end of inflation. From the above solution it can be 
shown that [27] 

-^ ~ (mt)2 < 1 . (36) 

Pkin 

Thus, when the vector field is light its energy density is dominated by its kinetic density. 
Therefore, [27] 



1 



Pa ^ Pkii 



^ (^end 



J^onrfA. 



end^^cnd 



'^end 



Pa oc a 



(37) 



i.e. the light vector field scales as radiation with the Universe expansion. 

When the physical vector field is heavy m ':^ H ^^ mt ^ 1 the solution in Eq. (34) 
becomes [27] 



A 



2 

77 



( l + 2d ^ 

Ci cos mt TT 



/ l-2d^ 

C2 COS mt TT 



(3^ 



which shows that the vector field is undergoing rapid quasi-harmonic oscillations whose 
envelope is decreasing as \\A\\ oc a~'^/^. This is easy to understand since, for a heavy vector 
field, within a Hubble time one can ignore the friction term in Eq. (32) and write it as 
A + m?A ~ 0. From Eqs. (33) and (38) it is straightforward to find 



1 2 2- 

Pa = —m t i+« 

IT 



cl + cl + 2ciC2 cos((i7r) 



Pa oc a" 



(39) 



where we used a oc t3(i+») in a spatially flat FRW Universe. Thus, we see that the density of 

the heavy oscillating vector field scales as pressureless matter with the Universe expansion. 

But is it pressureless indeed? From Eqs. (27), (33) and (38) we readily obtain [27] 



1 2 2_ 

p± = m. t i+« 

vr 



cl sin(2mt — dn) + Cg sin(2mt + dn) + 2ciC2 sin(2mt) =^ p± = , (40) 



i.e. over a Hubble time (which corresponds to a large number of oscillations) the average 
transverse pressure is zero. Since the longitudinal pressure is —p± it is zero too. This 



means that the energy-momentum of the rapidly oscillating homogeneous vector field is 
that or pressureless isotropic matter (cf. Eq. (26)). Hence, the vector field can dominate 
the Universe without introducing excessive anisotropic stress. One way of understanding 
this is that, due to the harmonic oscillations which send A — )■ —A, the direction of the 
vector field is rapidly alternated, so that, over a Hubble time, there is no net direction and 
the vector field behaves as an approximately isotropic fiuid. 



4.2 The vector curvaton scenario 

We have seen that the density of a massive Abelian vector field homogenised by infiation 
scales as radiation when the field is light and as matter when it becomes heavy. After the 
end of infiation the energy density is eventually transferred into a newly formed thermal 
bath of relativistic particles. The density of this thermal bath is dominated by radiation 
p^. The homogeneous vector field is initially light so its density also scales as radiation. 
Therefore, the density parameter VLa = Pa/ p of our vector field remains constant, where 
p ~ p^ is the density of the Universe. During radiation domination the Hubble parameter 
reduces with time as H{t) = l/2t so that eventually the vector field becomes heavy and 
begins its coherent oscillations. From then on, its density scales like matter and its den- 
sity parameter grows Qa oc a. Thus, the vector field has a chance to dominate (or nearly 
dominate) the Universe before its decay. When it does so it imposes its own curvature per- 
turbation onto the Universe, according to the curvaton scenario [19], without introducing 
any anisotropic stress [26]. A schematic representation of the vector curvaton scenario is 
presented in Fig. 1. 

The superhorizon perturbations of the vector field satisfy the same equation of mo- 
tion as Eq. (32). The reason is that this equation is linear and that the gradient term 
is heavily diluted for superhorizon perturbations.^^ Hence, the perturbations follow the 
same behaviour as the homogeneous zero- mode (with k = 0). Thus, when the vector field 
becomes heavy they undergo quasi-harmonic oscillations too and their anisotropic stress is 
also eliminated. 

The existence of these perturbations of the vector field implies that the density pa is 
also perturbed and the field's (near) domination of the Universe occurs at slightly different 
times at different locations. This results in a difference (perturbation) in the timescale of 
the Universe history, which, according to the SN philosophy, results in a contribution to 
the curvature perturbation (. Since the density pa is a scalar quantity, this is a scalar 
contribution to ( (and not a vector contribution). 

Let us now quantify the above. The curvature perturbation in the Universe is, in 
principle, the sum of the contribution of the vector field (a and any preexisting curvature 



^^In momentum space the gradient term is V^A — >■ {k/a)'^A where k/a is the physical momentum scale, 
which is k/a <C H for superhorizon perturbations. Note that the mass of the vector field is much larger 
than H, when it is oscillating and the same is true for its perturbations. 
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HOT BIG BANG 



Figure 1 : Log- log plot depicting the evolution of the density of the vector curvaton pA and the background 
radiation due to inflationary reheating p-y , from the end of inflation and until the onset of the hot big bang. 
During inflation, the vector field has a negligible contribution to the density of the Universe pA ^ Pinf • 
At the end of inflation (denoted by 'end'), the inflationary energy is given to a thermal bath of radiation 
Pini ~^ P-y (prompt reheating is assumed for simplicity). After the end of inflation p^ ex a~^, which is 
mimicked by the vector field density while the vector field remains light, i.e. pA oc a"**. Thus, the vector 
field density parameter V,a = Pa/ P remains constant with D.a ^ 1, where p = p-y + p^ is the density of the 
Universe (p ~ p^ during this period). At some later time (denoted 'osc') the vector field becomes heavy 
and begins oscillating. From then on, it behaves as a pressureless isotropic fluid whose density scales as 
PA oc a~^ . Thus, its density parameter grows as VLa oc a. This allows the oscillating vector curvaton to 
dominate the Universe at some later moment (denoted by 'dom'), when VLa — 1. Afterwards, the vector 
curvaton decays into the thermal bath of the hot big bang (its decay is denoted by 'dec'). The slanted 
dashed line corresponds to the possibility when the vector curvaton decays before domination, so that 
il'^'^ < 1. In this case the hot big bang begins at the original inflationary reheating. 
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perturbation, already present in the radiation fluid (^. Then we can write [19] 

c = C7 + a = (1 - ^aXj + ^aCa , (41) 

where CIa = ^_q — ^a and we have assumed that the vector field has already become 

heavy. In the above, Q corresponds to the curvature perturbation attributed to the i-th 
component of the Universe content, which, on a spatially fiat slice of spacetime, is given 
by |19] 

0.-ff^=l^, (42) 

Pi 3 Pi + Pi 

where we used the continuity equation pi + 3H{pi + Pi) = for independent fiuids. The 
above suggest that the contribution of the heavy vector field to ( is 

U = \aJ-<^. (43) 

3 Pa 

At the onset of the oscillations pkm ~ V^ so that p^ = 2V^ = ?Ti^||Ap. Thus, to first order 
we find^^ 

From Eq. (3) we see that the contribution of the vector field to ( to first order is (a = N\6Ai. 
Comparing with the above we get 

''A = pA^,^N^I = l^A%- (45) 

Thus, to second order, the contribution of the vector curvaton to the curvature perturbation 

is [13] 

2^ A6Ai , 1^ SAjSAj 

Note, however, that for 6 Ax/ A <^ 1, the one-loop correction (last term in the above) is 
negligible. 

Let us now turn our attention to non-Gaussianity. With the above values of N\ and 
N^ and Eqs. (22) and (23) it can be shown that [15] 



6 ,eql _ ,2^2 3 (1 + -2in + P+lip'- 2g^ 



A2 



^fZ = l3'ri-^^ '^\'' , ;-,, ^'^ (47) 



and 



6.sqz_«2^2 3 l + pAj + ipqA^yJl-Alsmu 



^The zero-mode and the perturbations begin oscillating simultaneously and in phase. 
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where p and q were defined in Eq. (21), A± is the projection of the unit vector A = A/A 
onto the plane of the three momentum vectors fci, ^2 and k^ which are used to define 
the bispectrum (cf. Eq. (14)) and u is the angle between A± and the shorter momentum 
vector in the squeezed configuration. 

From Eqs. (47) and (48) it is evident that, in the vector curvaton scenario with parity 
invariance, one can write 

/nl = f^l (l + GAl) (49) 

where Q is the anisotropy parameter for non-Gaussianity. In analogy to g in Eq. (2), Q 
quantifies statistical anisotropy in the bispectrum oi (. If ^ ^ 1, then non-Gaussianity is 
predominantly anisotropic, which means that /nl should have a clear angular modulation 
on the sky. If non-Gaussianity is indeed observed and no such modulation is found then 
models which predict Q ^ 1 will be ruled out. 

It is important to note above that the directions of statistical anisotropy in the spectrum 
and the bispectrum are correlated, since they are both determined by A. This is a smoking 
gun for the contribution of vector fields to (. 

5 Particle production of vector fields 

The vector curvaton mechanism can affect (or even generate) the curvature perturbation 
provided the vector field has, somehow, obtained a superhorizon spectrum of perturbations 
during inflation. In order to do this we need to have a mechanism which breaks the 
conformal invariance of the vector field. 

A massless Abelian vector field is conformally invariant, which means that it is not 
affected by the Universe expansion (it perceives it as a conformal transformation to which 
it is insensitive). Hence, it does not undergo particle production during inflation. Con- 
sequently, its quantum fluctuations do not give rise to classical perturbations of the field 
as is done, for example, with minimally coupled massless scalar fields. Thus, one expects 
a light vector field to be approximately conformally invariant and the production of its 
perturbations to be suppressed. An explicit breakdown of the vector field conformality is, 
therefore, required. This is model dependent, which suggests that observations might be 
able to discern between models and provide insight on the underlying theory. 

In this section we discuss two specific models which break the conformality of an Abelian 
vector field and have attracted considerable attention to date. Before going into these 
models however, let us discuss how we can use any such mechanism to obtain the spectra 
of perturbations of the vector field components. 

Assume, for the moment, that we are indeed operating under a suitable mechanism 
that breaks the conformality of the vector field. The first step is to perturb the vector field 
around the homogeneous value as A^ — > An{t) + SA^^x, t). Then we Fourier transform the 
perturbations as 

SA{k,t) = 1 5A{x,t)e~'^-''d^x 
and find the equations of motion of the Fourier components 5A{k^t) for the given model. 
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The next step is to promote the vector field perturbations to quantum operators by 
expanding in terms of creation and annihilation operators 

6A{x, t)= I -^ Y. [e-^cix{k)5A{k, t)e''' ' ^ + e^*d{{k)6A*{k, t)e~'^ ' ^] . (50) 

The mode functions 6A{k, t) and the Fourier components of the perturbations of the vector 
field satisfy the same equations of motion because the latter are linear. Thus, we need to 
solve these equations and find the the mode functions. To do this we need to employ the 
following boundary conditions, which are simply due to the fact that the perturbations 
begin as quantum fiuctuations well within the horizon [k/aH — )■ oo) 

pik/aH pik/aH 

5Al,r = -^=- and 6A\\ = 7 -^=^, (51) 

k/aH-*oc y/ 2k k/aH^oo y iK 



where 7 = — = \ [—Y + 1 is the Lorentz boost factor which takes us from the frame with 
A; = (where all components of the vector field perturbation are equivalent) to the one of 
momentum h. Apart from 7, we see that the vacuum boundary conditions are identical to 
the Bunch-Davis vacuum also employed for the particle production of scalar fields. 

Once we solve the equations of motion and find the mode functions we can obtain the 
power spectra of the superhorizon {k/aH -^ 0) perturbations using 

Vx = ^I^AP . (52) 

The typical value of the vector field perturbation is 6A\ ~ \/V\. Now, let us employ this 
method on two concrete models for the generation of a perturbation spectrum for the 
vector field during inflation. 

5.1 Non-minimal coupling to gravity 

This mechanism was first considered in Ref. [28] for the generation of a primordial magnetic 
field of superhorizon coherence. It was employed as a vector curvaton in Ref. [29] and also 
in Ref. [13]. 

Consider a massive Abelian vector field with a non-minimal coupling to gravity as 
follows 

C = -^F^uF'"' + ]pr? + aR)W^W^, (53) 

where R is the scalar curvature and a is a constant. The non-minimal coupling corresponds 
to a contribution to the effective mass of the field such that 

m^fj = TV? + aR . (54) 

The equations of motion have been found and solved for the above theory during de 
Sitter infiation, providing the following exact solutions for the mode functions 6A\ of the 
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perturbations. For the transverse components the solution is [29] 
5A.R = a-^'\ 



vr e'^^"-^^ 



Hl-e 



i2TTU 




(55) 



where 




- m ■ <-' 

The above solution produces a scale invariant spectrum if z/ = 3/2. This can be achieved 
if m -C if and a ~ |, because during de Sitter inflation R = —12H^}^ With this choice 
for m and a the solution for the longitudinal mode function is [13] 






^-(f)-' 



pik/aH 



2k 



(57) 



The transverse solutions are the same because the theory is parity invariant. However, 
it is clear that there is striking difference between the transverse and the longitudinal 
solutions. Using Eq. (52), we can now find the power spectra for the components of the 
perturbation. We obtain 

i.e. p = 1 and g = as expected (cf. Eq. (21)). Thus we find that particle production is 
anisotropic at a level of 100%. This means that the vector field contribution to ( should 
be subdominant, for otherwise it would violate the observational constrains which do not 
allow statistical anisotropy above 30%. Hence, we have to assume that ( is primarily due 
to some other source, presumably a light scalar field, and the contribution of the vector 
field is significant only at the level of generating significant statistical anisotropy. 
From Eq. (12), we obtain for the anisotropy parameter 

^ - /3<1 (59) 



1 + /3 



where we considered that V,f, = {H/2tiY and also that it is the scalar field which primarily 
modulates N so that /3 <^ 1 in Eq. (13). 

Using that p = 1 and g = in Eqs (47) and (48), we obtain 

Thus, we see that ^ ~ 1. This is because, in this theory, there is only one mass-scale 
involved, that is H. Dimensionless quantities, therefore, such as ^ or p are expected to be 



^■^This theory is in effect a modified gravity theory but it can be shown that the Friedman Equation is 
not affected if a = i and also RW^ is neghgible compared to the Einstein-Hilbcrt action for Wf^W^ <C ?7ip. 
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of order unity. In Ref. [15] it was shown that, whenever it is so, there is a clear prediction 
for the maximum non-Gaussianity, which provides a direct hnk with statistical anisotropy 
in the spectrum: 

/Nr~io'(^)'^'- (61) 

From the above, it is evident that, through the vector curvaton mechanism, significant 
non-Gaussianity can be produced. 

This theory was criticised in that it may suffer from instabilities such as ghosts [30]. 
However, it is not clear whether this is indeed so. In Ref. [30], it was shown that the modes 
of the longitudinal perturbations are ghosts but only when subhorizon. Given that these 
modes are subhorizon only for a limited time and also in view of the fact that the energy 
density of inflation is much larger than p^-, one may wonder whether these ghosts manage 
to destabilise the vacuum. A discussion on this issue can be found in Ref. [31]. 

5.2 Varying kinetic function and mass 

Consider a massive Abelian vector boson field, with Lagrangian density 

C = -\fF,,F^' + \m^W,W^, (62) 

where f = f{t) is the kinetic function, which is approaching unity by the end of inflation 
so that, afterwards, the vector field is canonically normalised. The above theory does not 
suffer from instabilities (ghost free) [32, 33], which motivates the use of the Maxwell-type 
kinetic term even if the field is not a gauge boson. Note here that an Abelian massive 
vector field is renormalisable even if it is not a gauge boson [34] . 

The spatial components of the physical canonically normalised vector field in this case 
are [27] 

A = ^W/a. (63) 

The mass of the physical, canonically normalised vector field is 

HI 

M^-^. (64) 

The massless version of this theory has been extensively considered, firstly for infla- 
tionary particle production for the generation of a primordial magnetic field [35] and more 
recently for the mild anisotropisation of inflation, which gives rise to statistical anisotropy 
in the curvature perturbation through anisotropic particle production of the inflaton [10]. 
Here we investigate if this vector field can play the role of the vector curvaton. This was 
first studied in Ref. [25]. 

The equations of motion for the mode functions for this model were obtained in Ref. [25]. 
They read 
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5Wl,r = (65) 



and 




5W|| = , (66) 



where 5W\ = a5A\/ \/J, cf. Eq. (63). Using the above, particle production has been 
studied in Ref. [27]. It was found that, the transverse components obtain a scale invariant 
superhorizon spectrum of perturbations when 

/ oc a~^^3 and M, < H^ (67) 

where the subscript '*' denotes the time of horizon exit. To obtain a scale-invariant 
spectrum for the longitudinal component one needs an additional condition on the time- 
dependence of ?n(t), which reads [27] 

m oc a . (68) 

If the vector field is a gauge boson then / is the gauge kinetic function which is related 
with the gauge coupling as / ~ 1/e^. This means that only the case when / oc a""^ is 
possible since, only then does the gauge field remain weakly coupled during inflation. The 
gauge kinetic function is one of the three fundamental functions which define a supergravity 
theory. ^^ In supergravity, the gauge kinetic function is a holomorphic function of the scalar 
fields of the theory. Now, during inflation, supergravity corrections are expected to give 
masses ~ if to the scalar fields [36]. This means that, during inflation, these scalar 
fields are fast-rolling down the slopes of the scalar potential, which would cause significant 
variation to the gauge kinetic function. Indeed, it is easy to show that /// ~ ii is natural 
to expect. Here we should note that, if / is modulated by the inflaton field then, under 
fairly general conditions, / oc a~^ is an attractor solution during inflation, which arises due 
to the backreaction of the vector field onto the roll of the inflaton down the inflationary 
scalar potential [37]. Thus, even though, originally, the inflaton may be also fast-rolling, 
the backreaction slows down its roll and allows slow-roll inflation to occur even with a 
relatively steep scalar potential. Indeed, it was shown in Ref. [37] that this is a neat way 
to overcome the infamous ?7-problem of inflation, while simultaneously obtaining / oc a~^ 
as an attractor solution. 

Now let us discuss the behaviour of the mode functions 5A\ in more detail. Firstly, we 

define 

k , M 

From the above and Eq. (64) it is evident that, if / oc a^ then z = constant, while if / oc a~'^ 
then z oc a^. Note that x oc a^^ always, while Eq. (67) requires that z* ^ 1. 

It turns out that there are three possible stages for the mode evolution [27]. When 
X > 1 ^ z, then the mode is still subhorizon and it is oscillating so that it can be matched 



The other two are the Kahler potential and the superpotentiaL 
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to the boundary conditions in Eq. (51). As time passes x decreases and the mode becomes 
superhorizon. When x, z <^ 1 then the mode is found to undergo power-law evolution. The 
third possible stage has to do with the case / oc a~^ only, when z is growing in time. In 
this case we could finally reach the time when z > 1 ^ x, when the superhorizon mode 
begins oscillating again. Since we are interested in superhorizon scales (they are the ones 
which can affect the curvature perturbation as observed in the CMB) we will consider the 
modes which are caught by the end of inflation either when x,z <^1, or when z > 1 ^ x, 
which is possible only with / oc a~^. 

5.2.1 Power-law regime 

This is the case when M < 3H <^ 2; < 1 when inflation ends. The mode functions for the 
superhorizon modes are found to be [27] 

SAlr = ^(^) and 6A\\ = -^(^)-. (70) 



Using this, the power spectra for the superhorizon perturbations of the vector fleld are 

where we used that the theory is parity invariant. Because M < 3H in this regime we flnd 
that P_|_ <^ V\\ , which means that particle production is strongly anisotropic. Therefore, 
if inflation ends in this power-law regime then the contribution of the vector curvaton to 
the curvature perturbation has to be subdominant (otherwise it would generate excessive 
statistical anisotropy), so that /3 ^ 1, which means that A^ is primarily modulated by a 
scalar fleld and not by our vector curvaton. In this case, therefore, the vector curvaton can 
only generate some statistical anisotropy in (. For the spectrum, the anisotropy parameter 
is (cf. Eq. (12)) 

9 = P^ = P/^', (72) 



where V^ = V+. Similarly, for the non-linearity parameter we obtain 

1 



t/nl — P 



5 2VL, 



1+ [p+^kp'^]A\ 



(73) 



where fi;=l{K = 0}for the equilateral {squeezed} conflgurations and we used Eqs. (47) 
and (48) and also that /3 <^ 1, g = and p = \j z^ (cf. Eq. (21)). Since the vector curvaton 
must have a subdominant contribution to C we have fi^ <^ 1, which gives ^a -^ ff^A- 
Using this and the above, we flnd that the isotropic part of non-Gaussianity has 

po _ 5^ _ 5^ .7. 
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while the anisotropy parameter for non-Gaussianity is 

^=p+ Kp2^1, (75) 

o 

i.e. non-Gaussianity is predominantly anisotropic. Thus, if non-Gaussianity is indeed 
observed and it does not feature a strong angular modulation this regime of this model 
will be ruled out.^^ 

If / oc a^ then M = constant and, because of Eq. (67), M = M^, <^ H. Thus, in this 
case Vj^,V\\ = constant and the above are the only possibihty. If, however, / oc a~^ then 
M (X. a^ and V\\ is decreasing in time. Yet, scale invariance is maintained because the 
amplitude of the modes when they exit the horizon is reduced in time in accordance to 
the reduction of the value of the spectrum. As a result, the piling of modes does not spoil 
the flatness of the superhorizon spectrum [27]. This case offers the possibility that M will 
reach and surpass H before the end of inflation. If this happens, then we are no more in 
the power-law regime. 

5.2.2 Oscillatory regime 

This regime corresponds to the possibility that M ^ H '^ z^ 1 when inflation ends. 
Because of Eq. (67), when the cosmological scales exit the horizon we have M* <C H. Thus, 
this regime can be realised only if M is growing during inflation, which is possible only if 
/ oc a~^. Note that, when M > H, particle production stops and the perturbation of the 
field on scales that leave the horizon becomes essentially zero (exponentially suppressed). 
However, the superhorizon scales which left the horizon when the vector field was still light, 
retain their perturbations, which evolve as follows. 
In this case, the mode functions are given by [27] 

L \ 3/2 • -1 

^ ^ '-^R = -^Jl/2{Z) + -x3v^J_i/2(^) (76) 

^yz 6 

-X^V^Jy2iz) - -^J_i/2{Z) . (77) 

Considering superhorizon modes (i.e. a; — )■ 0) and also that z ^ 1 in this regime, the above 
simplify to 

if\^/^sin^ , ^, 1 /i7\^/^cos2 

— and SAai = ^= — 




SAl,r = ^= [^] and 6A\\ = ^= — . (78) 

Thus, we see that the mode functions are undergoing rapid coherent oscillations since 
z ^ 1. Using these, we find that the average of the power spectra is 



""^^Note here that t/ ^ 1 because we have two scales into the model, H and M, in contrast to the 
non- minimally coupled model of the previous section, where we only had one scale H so that Q ^ I. 
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Thus, we see that we have isotropic particle production, which imphes p = q = (cf. 
Eq. (21)). In this case we do not need the input of a scalar field to generate the cur- 
vature perturbation as ( can indeed be fully produced by the vector field alone. To our 
knowledge this is the only model that manages to produce ( without the direct involvement 
of a fundamental scalar field. Because we need no scalar field, (3^1. Then Eqs. (47) and 
(48) both reduce to 

/nl = 4- , (80) 

which is identical to the case of a scalar curvaton [19]. 

5.2.3 Borderline regime 

It is interesting to briefiy consider the so-called "borderline" regime, when M ~ 3H -^ z ^ 1 
at the end of infiation. Again, this is realisable only if M is growing during infiation, i.e. 
when / oc a~^. 

In this case one has 

g = p=^-l = lL (81) 

where 6V = V\\ — V+. Thus, we see that, to satisfy observational bounds V\\ ~ V+ at least 
within 30%. If this is the case then it can be shown that [27] 

/nl = T^- (l + 9^1) . (82) 

Therefore, statistical anisotropy in the spectrum and the non-Gaussianity have the same 
magnitude, i.e. Q = g. This is an interesting characteristic signature for this scenario. 

5.2.4 The evolution of the zero mode 

When / oc a~^^^, it is straightforward to show that the equation of motion for the homo- 
geneous physical vector field is 

A + 3HA + APA = , (83) 

which looks identical to the Klein-Gordon equation for a minimally coupled massive scalar 
field. 

If / oc a^ then M = constant <^ 1 and we have A ^ constant. This means that 

Pa^Va-- MqAI = constant , (84) 

where Mq is the initial value of M (M = Mq, since M is constant) and Aq is the initial 
value oi A = \A\. 

Now, if / oc a~^ then M oc a^ and the solution to Eq. (83) is [27] 

A = Ao(^) V2cos (z±'^Y (85) 
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which means that the typical value of the vector field is A oc a ^. Using this, Eq. (33) gives 



Pkin 

which results in 



AqMq sin ( 2; ± 



AqMq cos Z ± 



and Va 



2 



^6) 



Pa = Pkin + Va = MqAI = constant . (87) 

Thus, we see that if / oc a~^^'^ then pA = constant during infiation. 

5.2.5 Curvaton physics 

Let us briefiy look into how the model parameters are constrained by the requirement 
that the vector field performs as a successful curvaton. As shown in Refs. [27, 32] such 
considerations impose the following constraint on the model parameters 

H^ " f~~, — f raaxiT A\ HdomW , n 

— ~ U^QT .\ "x , 88 

mp ^ \ mm{mA;-H*} / 



where 'dec' denotes the moment of the vector curvaton decay, 'dom' denotes the moment 
when the vector curvaton dominates the Universe (if it does not decay earlier than that), 
Ta is the decay rate of the vector curvaton and rriA is the final value of its mass at the end 
of infiation. 

If the vector field is still light at the end of infiation, i.e. if rriA <^ H^, we are in the 
anisotropic regime. Then we obtain [27] 

C~%^. (89) 

Using this and also the requirement that the curvaton decays before Big Bang Nucleosyn- 
thesis (BBN) (i.e. r^ > T^^^/mp where Tbbn ~ 1 MeV is the temperature at BBN) the 
constraint in Eq. (88) gives 

H^> y/^x 10^ GeV and 10 Tev <mA<^H^ . (90) 

If the vector field becomes heavy by the end of infiation, i.e. if niA ^ -f^*, then we are 
in the (almost) isotropic regime. In this case the vector curvaton alone can be responsible 
for C as we have discussed. Thus, we have [27] 

c = a ~ ^Tu => ^th. > enip . (91) 

This results in the following bound 

H, > 10^ GeV . (92) 

The above corresponds to a relatively high scale of infiation which, in supergravity models, 
may result in gravitino overproduction. However, the latter can be avoided through the 
entropy release by the vector curvaton decay which can dilute the gravitinos. 
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Taking into account that the vector curvaton can decay at least through gravitational 
couplings, we find F^ > m\/m?p. Using this it can be shown that the oscillations of the 
massive vector field cannot commence earlier than iV™^^ ^ 4.4 e-folds of infiation. As a 
result, the range of tua for which this regime can be realised is found to be [27] 

1 < !^ < io6. (93) 

From Eqs. (90) and (93) we see that there is ample parameter space for both the possibilities 
(light or heavy vector curvaton) to be realised. 

6 Conclusions 

The high precision cosmological observations enable cosmologists to investigate beyond the 
"vanilla" predictions of infiation and thereby discriminate between infiation models. A new 
such observable is statistical anisotropy, which amounts to direction dependent patterns 
in the CMB (and possibly large scale structure too). Statistical anisotropy is within the 
reach of the forthcoming observations of the Planck satellite mission, which is expected 
to release its first cosmological data in the beginning of 2013. Currently, the latest CMB 
observations allow statistical anisotropy in the spectrum as much as 30% {g < 0.3). Planck 
will reduce this bound down to 2% if statistical anisotropy is not observed. It should be 
noted here that, even if the spectrum is weakly statistically anisotropic, the bispectrum 
can be predominantly anisotropic with Q ^ 1. This means that, if non-Gaussianity is 
indeed observed without a strong angular modulation of /nl on the microwave sky, then 
all models which predict ^ S> 1 will be ruled out. 

Vector boson fields are natural candidates for the generation of statistical anisotropy 
in the curvature perturbation, because they are expected to pick a preferred direction 
when homogenised by infiation. The vector curvaton paradigm offers a simple, generic 
and effective mechanism for the direct contribution of vector boson fields to the curvature 
perturbation. The mechanism assumes a Proca vector field, whose zero-mode begins oscil- 
lating when the field becomes heavy, after the end of infiation. As shown, the oscillating 
zero-mode behaves as a pressureless isotropic fiuid and can (nearly) dominate the Universe 
without generating an excessive anisotropic stress. When doing so it imposes its own con- 
tribution to the curvature perturbation, in accordance to the curvaton mechanism. We 
should point out here that this perturbation is scalar in nature, because it is due to the 
perturbed value of the density of the vector field, which is a scalar quantity. A considerable 
advantage of the vector curvaton mechanism is that it does not rely on an interaction of 
any kind between the vector field and the infiaton sectors. This allows the vector field to 
correspond to physics at a much lower energy scale (e.g. TeV physics) than the scale of 
infiation, which may connect with observations in collider experiments such as the LHC. 

The particle production process, through which the vector curvaton obtains a super- 
horizon spectrum of perturbations, is in general anisotropic. This is because the vector 
boson field has more than one degrees of freedom (three if massive), for which the effi- 
ciency of the particle production is in general different. Thus, the curvature perturbation 
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contributed by a vector curvaton is, in general, statistically anisotropic. If particle pro- 
duction is indeed isotropic (at least at the level allowed by the observations) the vector 
curvaton mechanism can generate the curvature perturbation in the Universe from vector 
fields alone without directly involving any fundamental scalar fields. If, however, statisti- 
cal anisotropy is indeed observed, then we have to go beyond scalar fields to explain the 
observations. This means that, the observation of statistical anisotropy in the CMB, may 
probe the gauge field content of theories beyond the standard model. 

There are some related issues which we did not go into in this paper. For example 
studies of the trispectrum in vector field scenarios [14], or of one- loop contributions [12] 
Note also, the possibility that the vector curvaton is non-Abelian is investigated in Ref. [38] 
while the contribution to the curvature perturbation by P-forms can be found in Ref. [39] 

The interest now is in finding realistic candidates in theories beyond the standard model 
which can play the role of the vector curvaton. Examples can include the supermassive 
gauge bosons of grand unified theories [40]^^ or the vector fluxes on probe branes in the 
context of DBI- inflation [42]. Another promising idea is investigating the possibility that 
the vector bosons associated with gauged axions can play the role of the vector curvaton. 
In this case, the generation of parity- violating statistical anisotropy is possible [43]. 
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